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ABSTRACT 

o : 

S-^ ' We present an analysis of the stellar populations as a function of mass in a sam- 

c/3 ' pie of ^7000 galaxies of all morphological and emission types from the 6dF Galaxy 

Survey (6dFGS). We measure velocity dispersions and Lick indices from the spectra 
of the central regions of these galaxies, deriving ages and metallicities from the Lick 
, indices using stellar population models. We also derive dynamical masses and dynami- 

■ cal mass-to-light ratios for these galaxies by combining the central velocity dispersions 
(NJ ; with global photometry in the B, R and K bands from SuperCOSMOS and 2MASS. 
C^) . Together, these data allow us to reduce the degeneracies between age, metallicity and 

■ star formation burst-strength that have limited previous studies. 

We find that old galaxies exhibit a mass-metallicity relation with logarithmic slope 
. d[Fe/H]/dlogM » 0.25, while young galaxies show slopes consistent with zero. When 

' we account for the effects of the mass-metallicity relation, we obtain a single, consistent 

OO ! relation between mass-to-light ratio and mass for old galaxies in all passbands, M/L oc 

■ M ' 15 . As we have accounted for stellar population effects, this remaining variation in 
' the mass-to-light with mass (the residual 'tilt' of the Fundamental Plane) must have 

• i— I , a dynamical origin. However, we demonstrate that any simple trend between mass- 

■ to-light-ratio and mass or luminosity is inconsistent with the observations, and that a 
| more complex relationship must exist. 

We find that the central regions of galaxies of all masses often exhibit young 
stellar populations. However it is only in the lowest-mass galaxies studied (~ 10 10 M Q ) 
that these populations are evident in the global photometry. In higher-mass galaxies, 
young central populations have decreasing influence on the global photometry, with 
there being no discernible impact in galaxies more massive than - 2 x 10 11 M Q . 
We conclude that the young stellar populations detected in spectroscopic studies are 
generally centrally concentrated, and that there is an upper limit on the mass of star- 
forming events in massive galaxies. These results have important ramifications for 
mass-to-light ratios estimated from photometric observations. 

Key words: galaxies: general, galaxies: stellar content, galaxies: kinematics and 
dynamics 



1 INTRODUCTION 

Current models of galaxy formation and assembly can be 
characterised by two extreme and competing views. One 
such view is represented by the early "monolithic collapse" 
model, whereby the spheroidal (bulge) components of galax- 



ies form by the early collapse of individual gas clouds (e.g. 
Larson 1974; Carlberg 1984). In these early collapse mod- 
els the star formation induced by the collapse is terminated 
when the source of gas is either consumed or expelled from 
the galaxy. The bulge then ages passively (albeit perhaps 
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subject to the later addition of a disk). Such models are 
largely motivated by the small scatter in the observed rela- 
tionships between colour and magnitude (e.g. Baum 1959; 
Bower, Lucy & Ellis 1992a, b) and velocity dispersion, sur- 
face brightness and size (the Fundamental Plane; Dressier 
et al. 1987) for early-type galaxies and bulges. To explain 
the tightness of these relations, early collapse models require 
that star formation in bulges terminates at early times (z>2; 
Bower, Lucy &i Ellis 1992a, b), so that bulge dominated (e.g. 
elliptical) galaxies should be uniformly old. 

This picture faces challenges on a number of fronts. 
For instance, the differing slopes found for the Fundamen- 
tal Plane in various bands (e.g. Bernardi et al. 2003b) are 
difficult to explain in the monolithic collapse scenario. Fur- 
thermore, spectroscopic studies of the stellar populations in 
the bulges of galaxies find central ages ranging from ~1— 15 
Gyr (e.g. Trager et al. 2000; Proctor & Sansom 2002; Ter- 
levich & Forbes 2002: Caldwell, Rose & Concannon 2003; 
Kauffmann et al. 2003a,b; Proctor et al. 2004b), indicating 
evidence for more recent star formation. 

However, by far the most important challenge to the 
early collapse picture is the fundamental differences in 
galaxy formation and assembly histories that are predicted 
when the growth of structure via gravitational instability is 
modelled for a dark matter-dominated universe, now charac- 
terised extremely well via precision measurements of the cos- 
mological parameters (e.g. the LambdaCDM model; Spergel 
et al. 2003). Such "hierarchical" models, which represent the 
other extreme view, predict galaxy-mass objects to be as- 
sembled by the successive merging of lower-mass objects via 
hierarchical merging (White & Rees 1978), usually over an 
extended period. 

Another important prediction of hierarchical merging 
models is that the merger rate of galaxies is a function of 
environment, i.e. on average, galaxies in the dense environs 
of clusters merge early, while those in the less dense 'field' 
merge later (Kauffmann 1996). In the hierarchical merging 
scenario the properties of galaxies are therefore subject to 
strong influences from their environment. Results from stel- 
lar population studies appear to confirm the prediction of 
hierarchical merging models in that 'field' galaxies possess 
slightly younger central ages than their cluster counterparts 
(e.g. Thomas et al., 2005; Bernardi et al., 2006; Smith et al., 
2006). 

In contrast to early collapse models, hierarchical merg- 
ing models predict galaxy bulges to possess a range of 
ages, with low-mass galaxies forming earlier, and on shorter 
timescales, than high-mass galaxies. At first glance, the wide 
range of ages observed in spectroscopic studies might also 
appear to confirm the hierarchical merging predictions, how- 
ever the observed trend of age with mass is for lower-mass 
systems to be younger than higher-mass systems. This is in- 
consistent with hierarchical merging models which predict 
lower-mass systems to form earlier than high-mass system 
(Kauffmann 1996). The observed trend is therefore some- 
times referred to as 'anti-hierarchical'. 

The recent De Lucia et al. (2006) analysis of the Mil- 
lennium Simulation provides insights important to the de- 
velopment of a coherent picture of galaxy evolution. The 
Millennium Simulation models the evolution of dark mat- 
ter using an N-body code to which De Lucia et al. couple 
a semi-analytic model of galaxy formation. De Lucia et al's 



analysis clearly demonstrates that a strong distinction must 
be made between the mass assembly of galaxies and the for- 
mation of their stellar populations. This can be understood 
intuitively by the simple realisation that mergers bring old 
stars into a merger remnant as well as making new stars. 
Merger remnants may therefore be, and often are, still dom- 
inated (in terms of their mass) by old stellar populations. 
De Lucia et al. conclude that while the mass assembly of 
galaxies is indeed hierarchical, their star formation histories 
are nevertheless anti-hierarchical. Despite this improvement 
in our theoretical understanding, the problem remains how 
to reconcile the broad range of predicted (and spectroscopi- 
cally observed) galaxy ages with the extremely small scatter 
in some scaling relations. 

Recent surveys of large numbers of nearby galaxies 
have provided statistically significant results and hence fur- 
ther observational constraints on galaxy formation models. 
Such surveys include the 2dF Galaxy Redshift Survey (2dF- 
GRS: Colless et al. 2001b) and the Sloan Digital Sky Sur- 
vey (SDSS; York et al. 2000). From central region spectra 
and global optical colours the 2dFGRS probed the cosmic 
star formation history at relatively low redshift (Baldry et 
al. 2002) and its variation with local environmental density 
(Lewis et al. 2002). The SDSS also obtained optical colours 
and central spectra for a large number of low redshift galax- 
ies. This ongoing dataset has been used by several authors 
to study the star formation history and scaling relations 
of nearby galaxies. These studies include Bernardi et al. 
(2003b) who investigated scaling relations for ~9,000 galax- 
ies, Eisenstein et al. (2003) who combined spectra of 22,000 
massive galaxies, Gallazzi et al. (2005) who probed the star 
formation history of ~40,000 galaxies, Cid Fernandes et al. 
(2005) who conducted a spectral synthesis of ~50,000 galax- 
ies, and Chang et al. (2006) who combined spectra and op- 
tical colours with 2MASS near-infrared colours for ~3,000 
galaxies. Such large datasets have allowed these authors to 
investigate trends as a function of mass, environment, galaxy 
type etc. Already, a galaxy mass-metallicity relation has 
been detected in the metallicities of both gas (Tremonti et 
al. 2004) and stars (Gallazzi et al. 2005) and strong trends 
with galaxy environment identified (Lewis et al. 2002). 

However, the above studies have not explicitly explored 
the link between the central stellar populations, which they 
observe, and the global photometry of the galaxies in their 
samples. It is this issue that this paper specifically addresses. 
To this end, we compare the results of stellar population 
analysis of spectra from the 6dF Galaxy Survey (6dFGS; 
Jones et al. 2004) to the global near-infrared photometry of 
the 2MASS Extended Source Catalog (2MASS XSC; Jarrett 
et al. 2000a,b) and optical SuperCOSMOS (Hambly, Irwin 
& MacGillivray. 2001) data for several thousand galaxies. 
We explore such issues as the fraction of galaxies' stellar 
populations that can have been formed in recent star for- 
mation events and the effects of stellar populations on the 
slopes of colour- magnitude relations and the varying 'tilt' of 
the Fundamental Plane. 

The paper is laid out as follows. In the next section we 
outline the photometric and spectroscopic data. In Section 
[3] we detail the spectroscopic data reductions and analy- 
sis used to derive stellar population parameters. Section 2] 
presents the results of our analysis, which are discussed and 
summarised in Section [5] 
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2 THE SAMPLE DATA 
2.1 The 6dF Galaxy Survey 

The 6dFGS is a spectroscopic survey of the entire south- 
ern sky more than 10 degrees from the Galactic Plane. 
The survey measured redshifts for more than 124,000 
galaxies. The primary survey targets were selected pri- 
marily from the 2MASS Extended Source Catalog (XSC). 
All galaxies brighter than K toi =12.75 were included in 
the sample. Secondary samples of 2MASS and Super- 
COSMOS galaxies, complete down to magnitude limits 
(H, J, r F , bj) = (13.00, 13.75, 15.60, 16.75), were also in- 
cluded. 

The survey was carried out on the UK Schmidt Tele- 
scope using the Six-Degree Field (6dF) multi-fibre spectro- 
graph. The 6dF instrument simultaneously observed up to 
150 spectra over a 5.7° field. The science fibres of the instru- 
ment were 6.7 arcsec in diameter. The spectrograph is bench 
mounted in an enclosure inside the telescope dome, which 
gives it increased instrumental stability compared to a spec- 
trographs mounted on a telescope (e.g. 2dF and the SDSS 
spectrograph). This increased instrumental stability main- 
tains consistent instrumental resolution, which is important 
when working with Lick indices. The survey obtained spec- 
tra in the visual (V) and red (R) wavelength ranges provid- 
ing coverage from ~4000 to ~8400A. The spectra have 5-6A 
FWHM resolution in V and 9-12A resolution in R (see Jones 
et al. 2004 for more details). 

The main science aims of the survey included: to mea- 
sure the luminosity function of near-infrared selected galax- 
ies as a function of environment and galaxy type (Jones et 
al. 2006), to measure the clustering of galaxies, to produce 
detailed maps of the density and peculiar velocity fields of 
the nearby universe (Jones et al. 2004; Erdogdu et al. 2006) 
and to study the properties of the underlying stellar popu- 
lation of the galaxies such as ages and chemical abundances. 
It is this latter issue to which this paper is addressed. 

The 6dFGS is particularly suited to the aims of the cur- 
rent study because the selection criteria using 2MASS near- 
infrared photometry provides a sample less biased toward 
bright, young galaxies than surveys based on optical cata- 
logues. The K band photometry also has the useful property 
of being relatively insensitive to both metallicity the pres- 
ence of dust (unlike optical bands) , and is therefore an good 
tracer of mass. It also provides a useful tool in breaking 
the age-met allicity degeneracy. Finally, 6dFGS spectra are 
obtained from region 2 to 3 times larger on the sky than 
those used in the 2dFGRS and SDSS surveys. As a result, 
at the mean redshifts of the surveys, 2dF, SDSS and 6dFGS 
fibres correspond to 2.6, 3.9 and 4.6 kpc/h respectively. The 
6dFGS spectra therefore sample a significantly larger frac- 
tion of the total galaxy light than either 2dF or SDSS spec- 
tra. 



2.2 Spectroscopic data 

The spectroscopic data presented here are from the 6dFGS 
First Data Release (Jones et al. 2004). Lick index, veloc- 
ity dispersion and emission line-strength measurements 
were made on all galaxies with redshift quality > 3 
(the 'reliable' redshifts; Jones et al. 2004). This gave a 



total of 39649 galaxies with 3859 of these observed more 
than once (420 galaxies observed more than twice). If a 
galaxy was observed more than once, only the highest 
signal-to-noise-ratio spectrum was used for the measure- 
ments. However, repeat observations were used to establish 
the errors (see Section 13. 5p . Only galaxies with median 
signal-to-noise ratios >12 A -1 were deemed suitable for 
age and metallicity determinations. Such a value lies at 
the lower limit of signal-to-noise-ratio that can be used 
to derive ages and metallicities with precision. The above 
selection criteria resulted in ~7000 galaxies suitable for 
Lick index measurements. It should be noted that, as no 
morphological selection was made, the results presented 
here include galaxies of the full range of early and late types. 



2.3 Photometric data 

The K band magnitudes used in this work were taken to be 
the Ktot values estimated from the 2MASS data by Jones 
et al. (2004). We also use the re-calibrated SuperCOSMOS 
B and R band magnitudes (Hambly et al. 2001). Galaxy 
radii were taken as the 2MASS 20th magnitude per arcsec 2 
isophotal radius (Rk20\ Jarrett et al. 2001a). We also use 
K band half-light (effective) radii (R e //) for a spectroscop- 
ically selected early-type galaxy subsample, as supplied by 
Lachlan Campbell (private communication). Absolute mag- 
nitudes and physical galaxy radii were calculated using dis- 
tances based on the measured redshifts and a Hubble con- 
stant of 70 km s^ 1 Mpc -1 (our median redshift is 0.035). 
Errors on the distance estimates were assumed to be a com- 
bination of recession velocity uncertainty and peculiar mo- 
tions. In order to be conservative, we assumed the com- 
bined uncertainty to be ±1000 km s . The distributions 
of some key parameters in the final sample are shown in 
Fig. [1] The sample covers of order 5 magnitudes in lumi- 
nosity (-26 <Mjf< -21 mag and -22 <M S < -17 mag). 
This corresponds to a dynamical mass range of more than 
two orders of magnitude (from <10 10 to >10 12 Mq; Sec- 
tion \4A§. From the sample of ~7000 galaxies suitable for 
age-metallicity determination (Section 12. 20 ~6000 also had 
a full set of 2MASS and SuperCOSMOS photometry. 



3 SPECTROSCOPIC ANALYSIS 
3.1 Measuring emission lines 

We measured emission line-strengths for Ha, H/3, 
[NII]A6548, [NIIJA6584, [OIII]A4959 and [OIII]A5007 
from the 6dFGS spectra. In order to measure the Ha and 
[Nil] emission, a de-blending option in the IRAF task 
splot was used. This option allows the fitting of Gaussians 
to a series of spectral lines and the measurement of the 
value of the equivalent widths of each. The strengths of 
Ha and the two [Nil] doublet lines either side of it were 
measured in this way from continuum divided spectra. 
This process can lead to false measurements where noise 
spikes or the continuum level are fit with Gaussian profiles 
rather than real features. By consideration of the position 
and FWHM of the Gaussian fit to the feature, one can 
reasonably discriminate between real emission features and 
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at 4959 and 5007 A were measured using the procedure and 
band definitions of Gonzalez (1993). H/3 emission, on the 
other hand, was estimated for each galaxy as the difference 
between the observed 6dFGS H/3 index value and the value 
predicted by our best-fit model solutions (see Section l3~6l) . 

Within this work the sample is sometimes subdivided 
into 'passive' (no emission) and 'emission' galaxies. Passive 
galaxies are defined as those in which no emission lines are 
detected within the 6dFGS wavelength range with greater 
than approximately 1-cr significance. Emission line galaxies 
are therefore defined as those in which at least one emission 
line is detected at greater than one sigma confidence. For Ha 
this criterion corresponds to emission-line equivalent widths 
of ~1.5 A. Using this definition the sample contains 4640 
passive galaxies and 2270 emission galaxies. Unfortunately, 
a full analysis of morphological types is not available for 
the sample. However, visual inspection of a sub-sample of 
passive galaxies suggests that ~80% are early types, while 
amongst the emission galaxies a similar fraction are late 
types. 
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Figure 1. The distributions of key parameters in the final sample 
of ~7000 galaxies presented here. The signal-to-noise-ratio per 
A is estimated as the median of the 6dFGS spectra. Apparent 
magnitudes are from 2MASS and SuperCOSMOS data. Absolute 
magnitudes were calculated using redshifts from the 6dFGS spec- 
tra and a Hubble constant of 70 km s _1 Mpc _1 . K band 20th 
magnitude isophotal radius is from 2MASS. Effective radii are 
also presented for a (~25%) subsample of our 6dFGS galaxies. 
Dashed lines in the radius plots represent the radius of the 6dF 
aperture (3.35 arcsec). 



false detections. If the line centre was more than 7 A away 
from the features expected wavelength (given the 6dFGS 
redshift), the measurement was discarded. If the FWHM 
was well below the instrumental resolution of the 6dF 
instrument (< 5 A) the measurement was also discarded. 
Finally, extremely large FWHM (> 12 A) were discarded. 
These large FWHM are mostly Gaussian fits to features 
in the continuum, but may also have occasionally included 
real features of broad-lined active galactic nuclei (AGN). 
The equivalent widths of the two [OIII] emission lines 



3.2 Measuring velocity dispersion 

The method used to measure Lick indices from appropri- 
ately broadened spectra and to correct the indices for the 
broadening effect of a galaxy's velocity dispersion was based 
on that used by Stephen Moore (2001) 

Measurements of the velocity dispersion in the 6dFGS 
galaxy spectra were made using the IRAF task fxcor which 
carries out a Fourier cross-correlation between the object 
spectrum and a zero-redshift, zero- velocity dispersion tem- 
plate spectrum. The task outputs the radial velocity of the 
object spectrum and the FWHM of the Fourier correlation 
peak between the object and template spectra, which can 
be converted into a velocity dispersion. The template spec- 
tra used were three velocity standard stars observed with 
the 6dF spectrograph in 2003 March by Craig Harrison 
(HR2574, HR3145 and HR5888, which are late G and early 
K giants). 

The raw values of the velocity dispersion measured by 
fxcor need to be calibrated. To do this each of the stan- 
dard star spectra were artificially broadened to a specific 
velocity dispersion using the IRAF task gauss. These arti- 
ficially broadened spectra were then used to calibrate their 
(known) velocity dispersions against the value measured by 
fxcor. Towards lower velocity dispersions, as one approaches 
the instrumental velocity resolution (<100 km s _1 ), this cal- 
ibration becomes less accurate. However, galaxies with such 
low velocity dispersions require only very small corrections 
to their measured Lick index values. Consequently, reliable 
index determinations can still be obtained in low velocity 
dispersion galaxies. More importantly, this method of mea- 
suring velocity dispersion fails on galaxies that are not well 
fit by the standard star templates. This mainly occurs in low 
signal-to-noise ratio galaxy spectra with young ages and/or 
strong emission lines. 

A comparison of velocity dispersion values derived from 
the 6dFGS spectra to values in the literature (Hyperleda) 
is shown in Fig. [2] The 6dFGS values show good agree- 
ment with the literature values with an average offset of - 
6.1 kms -1 , and an RMS scatter of 20 kms -1 . Both deviation 
and scatter are greatest at high velocity dispersion, where 
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Figure 2. A comparison of velocity dispersion derived from 
6dFGS spectra (in km s — 1 ) with values from the literature (Hy- 
perleda). Despite the large 6dF aperture we find an offset from 
the one-to-one line (shown as a solid line) of only 6 km s with 
an RMS scatter of 20 km s _1 . Such differences have only small 
effects on Lick index determinations. 



we find the Hyperlead values to be consistently greater than 
the 6dFGS values. This is in agreement with Bernardi (2007) 
who found Hyperleda values greater than values derived 
from the SDSS at high velocity dispersion. 

Given the large 6dF aperture (6.7 arcsec compared to a 
more typical 2 or 3 arcsec in Hyperleda) it is possible that 
these deviations are the results of aperture effects. j0rgensen 
et al. (1995) investigated the effects of varying the aperture 
size in observations of velocity dispersions in the inner re- 
gions (i.e. R < H e ff/ 2 ) of early-type galaxies. Their results 
showed a weak power law dependence of measured velocity 
dispersions with the size of aperture. We therefore exper- 
imented with correcting velocity dispersions according to 
j0rgensen et al. (1995) in the sub-sample of 6dFGS galaxies 
for which Hyperlead values were also available. As this sub- 
sample of the data is heavily biased to nearby, early-type 
galaxies, it is well matched to the galaxies used to calibrate 
the J0rgensen et al. (1995) relation. The results showed that 
while the offset was reduced to +0.75 kms -1 , the scatter was 
increased to 23 kms -1 . We also note that the trend for Hy- 
perleda values to be greater than values derived from our 
data remains unaffected by these corrections. 

As noted above, the j0rgensen et al. (1995) correction is 
calibrated for the inner regions of early-type galaxies. How- 
ever, for our 6dFGS sample as a whole, the average aperture 
is ~1 R e // (see Fig.[T]). The sample also contains both early- 
and late- type galaxies. This suggests that the corrections 
presented by J0rgensen et al. (1995) are inappropriate for 
the 6dF sample as a whole. For such a sample, it is instead 
informative to consider the study of Gebhardt et al. (2000) . 
This analysis of the variation in observed velocity disper- 
sions (including the effects of rotation) measured through 
circular apertures of varying radius in a sample of 26 galax- 
ies of mixed Hubble types shows that, within 5 R e //, the 
sample exhibits little or no systematic variation of the mea- 
sured velocity dispersion with aperture size. 

As a result of the considerations above, we have elected 
to make no corrections for aperture size to the measured 
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velocity dispersion values within this work. 

3.3 Measuring Lick indices 

To make Lick index measurements suitable for direct com- 
parison with single stellar population (SSP) models, it is 
necessary to broaden the spectra to the instrumental reso- 
lution of the Lick/IDS system. The Lick/IDS FWHM instru- 
mental resolution as a function of wavelength is described in 
Worthey & Ottaviani (1997) where it is shown to vary from 
~8.5 A at the central wavelength (~5000 A) to ~10 A at 
the spectrum ends (~4000 and 6000 A). 

The calculation of the necessary broadening requires the 
measurement of the 6dF instrumental resolution. This was 
measured from 6dF arc line spectra obtained from the ob- 
servational time-span of the 6dFGS First Data Release. The 
instrumental resolution of the 6dF spectra varies not only 
with wavelength but also with optical fibre number, so mea- 
surements were made of arc lines at various wavelengths, 
for all fibres individually. A bivariate polynomial was fit to 
the measured FWHM spectral resolution as a function of 
both wavelength and fibre number. The amount of broad- 
ening required to match the Lick/IDS resolution could then 
be estimated and each 6dF spectrum was broadened to this 
resolution. 

3.4 Velocity dispersion correction 

The index values given in stellar population models do 
not include the effect of the velocity dispersion broaden- 
ing present in galaxy spectra. It is therefore necessary to 
introduce a correction to the observed galaxy index values 
to account for this effect. 

To measure the velocity dispersion correction required 
for each Lick index, we measured the indices of a range of 
artificially broadened stellar spectra. These values of the 
Lick indices measured in the broadened spectra were then 
compared to those measured in the un-broadened spectra 
to determine the index corrections for a range of velocity 
dispersions. 

For each index measured in Angstroms (e.g. the Fe and 
Ca indices) a multiplicative correction was defined by: 



Index correction = 



Index CT= 
Index 



(1) 



While, for each index measured in magnitudes (e.g. CNi, 
CN2, Mgi, Mg2), an additive corrections was defined by: 



Index correction = Index., 



Index™ 



(2) 



Where indexo- = o is the value of the index in the raw stellar 
spectrum (i.e. with no broadening) and index CT is the index 
value after broadening of the stellar spectrum to a specific 
velocity dispersion a. 

For each index, quadratic lines were fitted to the veloc- 
ity dispersion versus index correction ratio data. The index 
correction required for a particular galaxy could then be 
determined from the measured velocity dispersion and the 
quadratic function. 
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Figure 3. Top:The measured relationship between RMS error 
and signal-to-noise ratio from repeated observations for the G4300 
index (left hand plot) and for velocity dispersion (right plot). The 
small points are the individual differences between repeated ob- 
servations for each galaxy. The large points are the RMS binned 
values of these points after 2a clipping. The solid line is the least 
squares fit to the bin values. The dashed line is the relationship 
between error and signal-to-noise ratio after correcting for the 
clipping used in forming the RMS bins (see Section l3.5H The ver- 
tical line marks the lowest signal-to-noise ratio (12 A - 1 )used in 
the present study. Bottom: Same as top plots except the frac- 
tional error is shown. 



3.5 Estimating errors 

Estimating errors for spectroscopic measurements is notori- 
ously difficult, mainly because there are several sources of 
error that can not be well quantified (e.g. sky-subtraction, 
poor flux calibration etc). Therefore, instead of trying to 
propagate variance arrays, comparisons of the values from 
repeated 6dFGS observations of the same galaxy were used 
to quantify this error. This method was used for calculating 
the error in Lick indices and velocity dispersions, as well as 
in Ha and [Nil] emission equivalent-widths. 

In the First Data Release of the 6dFGS, 3859 galaxies 
were observed more than once. For all repeat observations 
the differences in values were grouped into bins based on 
the median signal-to-noise ratio of the galaxy spectra. The 
signal-to-noise ratio of the lower of the two observations was 
used for this binning, as this observation would be the cause 
of most of the error. If a galaxy had been observed more 
than twice then the two highest signal-to-noise ratio spectra 
were used. The RMS of each signal-to-noise ratio bin was 
determined and used as the o~ax value. To remove the biasing 
effect of large outliers in the bins only the lowest 95.45% of 
the values (the 2a limit) in each bin was used to calculated 
the RMS value. This requires that the RMS be scaled up by 
1.164 to take into account this clipping. 

A function was fitted to the RMS of the bins against the 
signal-to-noise ratio. The fitting function used was a modi- 
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Figure 4. Mean residuals of indices to the final best-fit SSP 
models are given in units of observational error (i.e. \). Error 
bars represent the RMS scatter about the mean. Solid symbols 
represent indices included in the fits, while open symbols repre- 
sent those indices excluded from the fits. The results for fits to 
both passive (left) and emission galaxies (right) are shown. The 
average reduced— \ 2 of indices included in each subsample (x?) i s 
indicated in each plot. 



lied rectangular hyperbola of the form y — a/(x + b) ('a' and 
'b' the fitted coefficients). The error in any measured quan- 
tity for a galaxy could then be estimated using the signal- 
to-noise ratio of the spectrum and the corresponding error 
from the fitted function. Examples of this procedure can be 
seen in Figure [3] We note that the majority of the sample 
with repeat measurements were re-sampled due to the low 
signal-to-noise of the original observation. 



3.6 Fitting indices to models 

We used the x 2- fitting procedure of Proctor & Sansom 
(2002) (see also Proctor et al. 2004a,b and Proctor et al. 
2005) to measure the derived parameters: log(age), [Fe/H], 
[Z/H] and [E/Fe] (a proxy for the 'a '-abundance ratio; see 
Thomas et al. 2003 for details). Briefly, the technique for 
deriving these parameters involves the simultaneous com- 
parison of as many observed indices as possible to models 
of single stellar populations (SSPs). The best fit is found 
by minimising the square of the deviations between obser- 
vations and models in terms of the observational errors (i.e. 
X 2 )- The rationale behind this approach is that, while all in- 
dices show some degeneracy with respect to each of the de- 
rived parameters, each index does contain some information 
regarding each parameter. In addition, such an approach 
should be relatively robust with respect to many problems 
that are commonly experienced in the measurement of spec- 
tral indices and their errors. These include poor flux calibra- 
tion, poor sky subtraction, poorly constrained velocity dis- 
persions, poor calibration to the Lick system and emission- 
line contamination. This robustness is of particular impor- 
tance in the analysis of large numbers of pipe-line reduced 
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spectra such as those of the 6dFGS which cannot be flux 
calibrated and so are not fully calibrated to the Lick sys- 
tem. The method is also relatively robust with respect to 
the uncertainties in the SSP models used in the interpre- 
tation of the measured indices; e.g. the second parameter 
effect in horizontal branch morphologies and the uncertain- 
ties associated with the Asymptotic Giant Branch. It was 
shown in Proctor et al. (2004a) and Proctor et al. (2005) 
that the results derived using the x 2 technique are, indeed, 
significantly more reliable than those based on only a few 
indices. 

In order to carry out the comparison of observations 
to SSP models it was first necessary to select and inter- 
polate the models provided in the literature. Interpolations 
are required because models in the literature are presented 
at only five or six discrete metallicities over the ~2.5 dex 
range. We experimented with a number of different models 
(Bruzual & Chariot 2003 (hereafter BC03); Thomas, Maras- 
ton & Bender 2003; Thomas, Maraston & Korn 2004 and 
Korn, Maraston & Thomas 2005 (hereafter KMT05). The 
results presented here are those from the models of KMT05, 
as these provide the necessary coverage in age (0.1 to 15 
Gyr) and metallicity ([Z/H] from -2.3 to +0.67), as well as 
including a detailed modelling of the effects on indices of 
varying 'a'-abundance ratios. However, fits were also ob- 
tained to the other model sets for comparison. 

The process by which the best fits were obtained was it- 
erative. First, fits were obtained for all galaxies using all the 
available indices. The pattern of deviations from the fit so 
obtained was then used to identify indices that matched the 
models poorly (see Fig. Due to the problem of emission- 
line filling of the H<5, H7 and H/3 Balmer lines in emis- 
sion galaxies, we considered 'passive' and 'emission' galax- 
ies (Section [37TJ separately. In Fig. Remission galaxies show 
the expected large deviation between observed and best fit 
model values of Balmer indices and Fe5015 (which is also 
emission- line affected). These indices were therefore omit- 
ted from the fitting procedure in emission galaxies. 

Perhaps surprisingly, H/3 and Fe5015 also deviate sig- 
nificantly in the 'passive' galaxies. We take this to indicate 
that, despite the rather severe definition of emission galaxies 
outlined in Section ^. II at least some of the passive galaxies 
contain low-level emission that has escaped detection during 
our emission line measurements. Indeed, this was the main 
motivation for the stringent definition of emission-galaxies. 

To minimise the effects of this low-level emission on our 
age and metallicity estimates, age sensitive H/3 and metal- 
licity sensitive Fe5015 indices were also excluded from the 
fits to passive galaxy data. 

In both passive and emission galaxies Ca4455 was found 
to fall below best-fit model values. This effect has been noted 
in previous studies (e.g. Proctor & Sansom 2002; Clemens 
et al 2006). The flux-calibration sensitive Mgi and Mg2 in- 
dices were also found to deviate from best fit model values. 
This was an anticipated effect, as 6dFGS spectra are only 
approximately flux calibrated and Mgi and Mg2 have widely 
space pseudo-continuum bands. Ca4455, Mgi and Mg2 were 
therefore excluded from the fitting procedure. Fortunately, 
most of the information lost by the exclusion of the Mgi 
and Mg2 indices is captured in the Mgb index that is not 
sensitive to flux-calibration issues. 

Once poorly fitting indices had been identified and 
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Figure 5. The [OIII]/H/3 vs [Nil] /Ho emission line diagnostic 
diagram. Nil, Ha and [OIII] were measured directly (see Sec- 
tion [37TJ. H/3 emission was estimated as the difference between 
observed and best-fit SSP model values of the H/3 index (see Sec- 
tion I3.6H . The average errors in the line-ratios are shown in the 
bottom-left of the plot. The solid line is from Kewley et al. (2001) 
and represents the theoretical dividing line between star-burst 
galaxies (grey points) and AGN (black points). The dashed line 
is the equivalent line from KaufTmann (2003a). Conservatively, 
we use the Kewley line to identify AGN. 



omitted, the fitting of the remaining indices was carried out 
using a clipping procedure in which indices deviating from 
the model fit by more than 3<r were excluded, and the fitting 
procedure performed again. This results in an average of less 
than one additional index per galaxy being clipped from the 
fits. Most of the clipped indices could be associated with 
contamination from features such as the 5577 A sky-line. 
The resultant fits are based on between 10 and 15 indices 
and show good agreement between observed and best- fit val- 
ues (Fig. |4}. Indeed, the reduced-^ 2 values given in Fig. [4] 
(Xr) were 1.08 and 1.30 for passive- and emission-galaxies 
respectively. This clearly indicates that our data are well 
represented by the models and that our errors reflect well 
the differences between observations and models. 



3.7 Errors in derived parameters 

For each galaxy in the sample, errors in log(age) and [Z/H] 
were estimated using 50 Monte-Carlo realisations of the 
best-fit model indices, assuming the observational errors for 
the galaxy in question. We again note that the reduced-x 2 
values of both passive- and emission-galaxies were close to 
one, indicating that the assumed errors reflect the average 
differences between observations and models quite well. The 
RMS scatter in the log(age) and [Z/H] values of the 50 real- 
isations was therefore taken as the individual galaxy error. 
However, in an effort to better understand the detailed be- 
haviour of the errors, deviations from the best-fit log(age) 
and [Z/H] values of all 50 realisations of each galaxy were 
retained and are presented in Section \4. 31 For the purposes 
of the following analysis it is important that the derived ages 
are well defined. Therefore, ~ 1000 low signal-to-noise galax- 
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ies whose error in log(age) exceeded 0.3 dex were deemed 
unreliable and were excluded from the analysis. The aver- 
age errors for the remaining galaxies were ~0.15 dex in both 
log(age) and [Z/H]. 



4 RESULTS 

4.1 Emission line galaxies 

Emission lines can be used to separate star-burst and active 
galactic nuclei (AGN) galaxies using the [OIII] A5007/H/3- 
[NII]/Ha diagnostic. Values of the [OIII]A5007, [Nil] and 
Ha emission strengths were estimated from the 6dFGS 
spectra in a traditional manner (Section 13. ip . The strength 
of H/3 emission, however, was taken to be the deviation 
between the measured value of the H/3 index and the 
best-fit SSP model value. Due to the relatively small 
observed dynamical range of H/3 in absorption compared to 
that observed in emission, this is a fairly robust estimate 
of H/3 emission. This is evidenced by the classic 'Y'-shape 
evident in the diagnostic plot of [OIII]A5007/H/3 against 
[Nil] /Ha shown in in Fig. [5] The position of the line 
dividing star-burst galaxies from AGN varies between 
studies in the literature. In Fig. [5] we show the lines given 
by Kewley et al. (2001) (solid line) and Kauffmann (2003a) 
(dashed line). In order to ensure reliable identifications of 
AGN we use the definition of Kewley et al. in the following. 
About 2200 of the 2700 emission galaxies in our sample 
have reliable estimates of all four of these emission lines, 
300 of these lying in the region of the diagram associated 
with AGN as defined by Kewley et al. (2001). 



4.2 Age and metallicity 

In this section we outline the results of our age and metal- 
licity determinations using Lick indices. 

The distribution of ages and metallicities are presented 
in Fig. [6] The sample has been sub-divided into passive- 
and emission-galaxies as described in Section T3, II Emission- 
galaxies classified as AGN (Section 14. ip are shown as grey 
symbols, the remainder are plotted in black. 

Before interpreting this plot, we must note the apparent 
clusterings and 'zones-of-avoidance', the clearest examples of 
which are evident in the emission galaxy plot. These are the 
inevitable result of linear interpolations in non-rectilinear 
spaces. 

Returning attention to the actual results of our age 
and metallicity determinations; similar trends of increasing 
[Z/H] with decreasing log(age) are evident in both of passive 
and emission galaxies - albeit with emission galaxies tend- 
ing to younger ages. Galaxies classified as possessing AGN 
in Section ^. ll also often possess young central ages. 

We investigate the effect of aperture size on derived pa- 
rameter using Fig. [7] In these plots, the galaxies with the 
largest and smallest sizes (Ra'20 >37.5 arcsec and Rk20 <7.5 
arcsec) are shown as red and blue points respectively. The 
blue points therefore represent the galaxies in which the 6.7 
arcsec AAOmega aperture encompasses the largest fraction 
of galaxy light (on average ~3.0 R e //), while the red points 
represent galaxies in which the smallest fraction of galaxy 



light (~0.5 R e //) is encompassed by the fibre. The observed 
age-met allicity and age-velocity dispersion relations can be 
seen to be similar in these two extremes, with the most no- 
ticeable difference a slightly higher [Z/H] (by ~0.1 dex) in 
large galaxies. We therefore find aperture-size effects to be 
generally small. 

Due to the large, circular aperture of the 6dFGS, a di- 
rect comparison with the literature, which are mainly based 
on long-slit observations of galaxy centres, is not straight 
forward. Nevertheless, we note that the observed trend is 
at least qualitatively in agreement with the trends observed 
in previous studies (Trager et al. 2000; Proctor & Sansom 
2002; Mehlert et al. 2003; Gallazzi et al. 2005; Collobert et 
al. 2006). The young central ages of AGN are also consis- 
tent with studies in the literature showing that star-bursts 
are often found to accompany AGN activity (e.g. see Cid 
Fernandes et al. 2004 and references therein). 

An intriguing feature of Fig. [5] is the presence of a 
few old galaxies with extremely low [Z/H] values, (<-0.5 
dex). These galaxies, present in both passive and emission 
samples, possess a range of velocity dispersions. Visual in- 
spection of both spectra and photometry generally reveals 
nothing exceptional about these galaxies. Further analysis 
of these galaxies are beyond the scope of this paper, but we 
note that their low numbers result in them having no impact 
on the results and conclusions of the present paper. 

As in these literature studies, the slope in the age- 
metallicity relation exhibited by our data is similar to the 
slope of the age-met allicity degeneracy (Worthey 1994; 'the 
3/2 rule'). Before comparing these results with the photom- 
etry, it is therefore clearly important to establish that the 
age-met allicity degeneracy has been broken. To this end, in 
the following section we consider the errors in the derived 
parameters. 



4.3 Errors in age and metallicity 

As described in Section 13.71 errors in age and metallic- 
ity were characterised by considering the results of the 50 
Monte-Carlo realisations of each of the 7000 galaxies in our 
sample. Fits using both the index combinations shown in 
Fig- E] (with and without Balmer lines) were carried out 
for the purposes of this analysis (a total of approximately 
one million realisations). The differences between the input- 
model values of log(age) and [Z/H] and those of each of the 
realisations were calculated. These were combined in a num- 
ber of bins depending upon the age-metallicity of the input 
models. 

The analysis was also carried out using the observed 
index values as the inputs to the Monte-Carlo realisations 
(rather than the best-fit model values). No significant quan- 
titative or qualitative differences were found. 

Fig. [8]shows the results of this analysis as 1-sigma confi- 
dence contours for three age-metallicity bins. Each contour 
in these plots is based on tens of thousands of individual 
realisations of the best-fit galaxy data. The contour lev- 
els were defined as the iso-densities corresponding to e _1 
of the peak values (as expected for the 1-sigma contour of 
a two-dimensional Gaussian distribution). The use of the 
iso-density peak in estimating the extent of 1-sigma con- 
fidence contours results in a small uncertainty in the esti- 
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Figure 6. The ages and metallicities derived from the 6dFGS spectra separated by emission type (see Section [37TJ . Grey points in the 
right-hand panel represent the AGN identified in Fig. [5] Average errors are shown in each plot. The overall loci of the points show 
that passive and emission galaxies largely follow similar trends between age and metallicity - albeit with emission galaxies tending to 
significantly younger ages. 



mates. However, we estimate these to be of order 2% and 
therefore insignificant. 

Results of the error analysis for both passive galaxies 
(in which Balmer indices were included in fits) and emission 
galaxies (in which Balmer indices were excluded from fits) 
are shown separately. We recall that H/3 was excluded from 
all the fits. Fig.[S]also shows the data divided into two signal- 
to-noise ratio regimes (S/N<16 A" 1 and S/NM6 A -1 ). 
Marginal distributions in log(age) and [Z/H] are shown at 
the edges of each main plot. We note that the positions of 
the contours in Fig. [8] were chosen for clarity only. To get 
a full picture one must imagine the whole surface of the 
log(age)-[Z/H] plane populated by such contours. 

Fig.|8]shows that, as might be expected, error estimates 
increase with decreasing signal-to-noise-ratio. For passive 
galaxies, the analysis reveals log(age) and [Z/H] errors of 
~0.1 dex in galaxies with S/N above ~16. This falls to 
~0.2 dex in galaxies with S/N below ~16. The errors of the 
emission-galaxies (in which all Balmer lines are excluded 
from the fits) are somewhat larger, but are generally still 
of order 0.2 dex. It is therefore clear that while the age- 
metallicity degeneracy is still present (as evidenced by the 
sloped elliptical error contours in Fig. [Sj , the magnitude of 
the errors is sufficiently small for its effects to be negligible 
in our results - even in galaxies fit without Balmer lines. 
This is the result, and main advantage, of using large num- 
bers (>10) indices in the determination of ages and metal- 
licities. The good age resolution achieved is emphasised in 
the marginal distributions of Fig. [8] in which galaxies older 
than 10 Gyr (1.0 dex), and galaxies younger than 1.5 Gyr 
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Figure 7. Age-metallicity and age-log(cr) are shown with colour 
representing galaxy (apparent) size. Blue points represent the 
smallest 5% of our sample, red points the largest 5%. Age- 
metallicity and age— log(cr) are shown to be largely independent 
of galaxy size. 



(0.15 dex) (delimited by dashed lines in the marginal distri- 
butions) are largely uncontaminated by galaxies with ages 
of 3 Gyr (0.5 dex). Since the following analysis concentrates 
mainly on the very oldest and very youngest galaxies, these 
provide confidence that results do not suffer significantly 
from the effects of the age-metallicity degeneracy. 

The final sample for which velocity dispersions and 
reliable age/metallicity estimates were measured consists of 
4500 passive galaxies and 1000 emission galaxies. 
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Figure 8. Error contours derived from the Monte-Carlo realisations of our galaxy sample are shown in three age-metallicity bins. Two 
signal-to-noise ratio subsamples are shown; S/N>16 A" 1 (left) S/N and <16 A" 1 (right). Results are given for both passive- (top) and 
emission-line (bottom) galaxies using the appropriate index combinations (see Section l3.6l and Fig.Q. The data for galaxies of appropriate 
emission characteristics and signal-to-noisc-ratio are shown as grey points in each plot. Dashed lines in marginal distributions represent 
the upper and lower limits used to define the 'young' and 'old' samples respectively. Although the age— metallicity degeneracy is still 
present, its effects are minimised by the relatively small errors. 



4.4 Photometric properties and dynamical masses 

The combination of 6dFGS, 2MASS and SuperCOSMOS 
data permits the comparison for spectroscopic age and 
metallicity determinations to the B, R and K band pho- 
tometric data for some 6000 galaxies. Specifically, by com- 
bining the photometry with the velocity dispersion measures 
from 6dFGS spectroscopy, we are able to investigate trends 
with age and metallicity in B, R and K band mass-to-light- 
ratios. 

This analysis is carried out using a dynamical mass 
(M^yn) calculated as: 



M dyn = 



Co R c ff 



(3) 



Where a is the central velocity dispersion, R g is the half- 



light radius and the constant C has a value of 5.0 (Cappellari 
et al. 2006). Unfortunately half-light radii are only available 
for 25% of our sample, while K band 20th magnitude isopho- 
tal radii (Rk2o) are available for the entire sample. A cali- 
bration was therefore carried out by using Equation[3]to cal- 
culate masses for the galaxies in the 25% sub-sample (which 
also possess much more accurate, aperture corrected, veloc- 
ity dispersion estimates), and comparing them with mass 
estimates using Rk20 (instead of R e ff ) in the same equation. 
A plot of the comparison is shown in Fig. [9] The correlation 
has equation Md yn =1.068MR K20 -1.223 with only 0.12 dex of 
scatter. Given error estimates in Md y n and Mr K20 of 0.06 
and 0.10 dex respectively, the correlation is clearly extremely 
good. Masses quoted throughout the remainder of this work 
are therefore those based on Rif20 corrected as detailed. 
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Table 1. Examples of B, R and K magnitudes from BC03 for 
single stellar population models of varying age and metallicity are 
presented. Note that these values are used almost exclusively in 
a differential manner (see Section 1 1. Il l . 



As well as the photometry of 2MASS and SuperCOS- 
MOS, we also use the photometric results from the SSP 
models of Bruzual & Chariot (2003). Some key values from 
the Bruzual & Chariot models are shown in Table [TJ We 
use these models almost exclusively to estimate differential 
properties; e.g. the rate of change of mass-to- light-ratio with 
log(age) in populations of given [Z/H]. It should be noted 
that the rates of change of the photometric properties are 
reasonably constant with metallicity and age, as long as age 
is expressed in logarithmic form. This clearly indicates that 
a differential approach is a robust use of the model values. 
We shall flag the one occasion in which the models are used 
in a non- differential manner. 



4.5 The effects of age and metallicity on 
photometric properties 

This section details how the ages and metallicities derived 
from 6dFGS spectra for the central regions of galaxies are 
related to their global photometry. The analysis considers 
mainly the youngest (age<1.5 Gyr) and oldest (age>10 Gyr) 
galaxies. By concentrating on the two extremes in age we 
minimise the effects of the age-burst-strength degeneracy 
on our age and metallicity estimates. However, it is also im- 
portant to note that the selection of galaxies with extremely 
young central populations has the effect of biasing the sam- 
ple towards the largest bursts of star formation. 

The following analysis utilises average values of galaxy 
parameters binned according to either mass or luminosity. 
The derived values are shown in Table [2] All colours and 
mass-to- light-ratios for the binned data were calculated from 
values in this table. We assume M.8,0^5.47, M_r,0=4.28 and 
Mx,©=3.33. 

A mass-metallicity plot of the two sub-samples of our 
data is shown in Fig. HUM . In this plot, galaxies younger than 
1.5 Gyr are shown as black symbols; these galaxies possess 
an average log(age) of 0.0 dex (1 Gyr). Galaxies older than 
10 Gyr are shown as grey symbols; these possess an average 
log(age) of 1.105 dex (12.7 Gyr). Results of the mass binning 
of metallicity values are shown in Fig. [10] as blue and red 
lines respectively. It is evident that the mass-metallicity re- 
lation detected varies with age. The logarithmic slope of the 
relation in old galaxies is ~0.25, while the young galaxies 
are more metal rich than the old, and the data are consis- 
tent with no slope. An important aspect of our analysis will 
therefore be a consideration of the effects of these trends 
on mass-to-light-ratios. We also investigate implications for 
optical/near infrared colour-magnitude relations. 



4.6 Mass-to-light-ratios in old galaxies 

We consider first the ~1500 galaxies identified by our spec- 
tral analysis to possess ages older than 10 Gyr. Plots of 
[M/L] against dynamical mass and K band luminosity for 
these galaxies are shown in Fig. 1111 In each plot, the average 
values of the data in five bins along the x-axis (Tabled are 
shown as solid lines. The extent of the rms scatter in each 
bin is identified by the dashed lines. 

The observed trends between mass and mass-to-light- 
ratio in old galaxies vary with both waveband and variable 
against which they are plotted. In the K band, we find no 
trend in mass-to-light-ratio with luminosity, i.e. a slope con- 
sistent with zero. The other five plots all show trends with 
varying, non-zero slope. The logarithmic slopes in mass-to- 
light-ratio with mass are 0.231, 0.214 and 0.157 in B, R and 
K bands respectively, while the logarithmic slopes with lu- 
minosity are 0.093, 0.080 and 0.022 respectively. The typical 
formal errors on these logarithmic slopes are ~0.01. These 
results are in good agreement with Trujillo, Burkert & Bell 



1 No significant difference was found between emission-types in 
these, and subsequent, figures. We therefore make no distinction 
between them, although we note that the majority of young galax- 
ies are also emission galaxies, while, conversely, the majority of 
old galaxies are passive galaxies (Fig. [51. 
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10.778 


-0.136 


1.105 


-19.308 


-20.395 


-23.530 


10.735 


-0.138 


1.107 


-19.027 


-20.112 


-23.297 


10.323 


-0.276 


1.101 


-18.529 


-19.546 


-22.623 


10.330 


-0.213 


1.104 


-18.232 


-19.222 


-22.340 


Young, Mass binned 










Young, Luminosity binned 








log(Mass) 


[Z/H] 


log(age) 


M fl 


Mh 


U K 


log(Mass) 


[Z/H] 


log(age) 


M s 


M R 


M K 


11.648 


0.363 


0.062 


-21.191 


-22.216 


-25.557 


11.430 


0.317 


0.042 


-21.366 


-22.382 


-25.728 


11.331 


0.421 


-0.033 


-20.745 


-21.733 


-24.959 


11.209 


0.311 


-0.017 


-20.781 


-21.711 


-24.948 


11.062 


0.394 


-0.054 


-20.432 


-21.317 


-24.569 


10.870 


0.276 


0.023 


-19.956 


-20.878 


-24.079 


10.781 


0.381 


-0.011 


-19.804 


-20.694 


-23.897 


10.563 


0.285 


-0.025 


-19.204 


-20.015 


-23.259 


10.353 


0.320 


-0.025 


-19.285 


-20.130 


-23.229 


10.170 


0.159 


0.028 


-18.545 


-19.303 


-22.332 



Table 2. Binned average values of mass, metallicity, age and luminosity of old (>10 Gyr) galaxies (top) and young (<1.5 Gyr) galaxies 
(bottom). Binning was carried out in both mass (left) and K band luminosity (right). 




10 11 12 



log(Mass) 

Figure 10. The relations between metallicity ([Z/H]) and mass 
is shown for young (<1.5 Gyr; black symbols) and old (>10 Gyr; 
grey symbols) galaxies. Mass-binned average values from Table 
[2]l are shown for both old galaxies (red solid line) and young 
galaxies (blue solid line). While old galaxies exhibit a clear mass- 
metallicity relation, young galaxies do not. 

(2004) who use a combined SDSS/2MASS catalogue to in- 
vestigate mass-to-light-ratios in the B and K bands (but 
without the luxury of spectroscopy). They found similar 
slopes in B and K band stellar mass-to-light-ratios with lu- 
minosity of 0.07 and 0.02 dex respectively. 

While the effects of age are almost entirely eliminated 
by our consideration of only old galaxies, the effects of 
metallicity have yet to be accounted for. The effects of 
the mass-metallicity relation shown in Fig. [10] on mass- 
to-light-ratios are investigated in Fig. [12] This figure again 
shows the binned data from Fig.[TT]as thick solid lines. Also 
shown are the values after correction of the mass-to-light- 
ratio for differences in metallicity (dashed lines). The cor- 
rections are made by comparing the luminosities predicted 
by BC03 models for SSPs of the metallicity (and age) of 
each mass or luminosity bin (Table [2]| with the correspond- 
ing prediction for the highest mass/luminosity bin (Table[3]). 
The dashed lines in the left hand plots therefore represent 
the observed trends between dynamical mass and mass-to- 




10 10.5 11 11.5 12 10 10.5 11 11.5 12 
log(M dyn ) logd/Lg,) 

Figure 11. The dynamical mass-to-light-ratios of old galaxies in 
B, R and K bands are plotted against dynamical mass and K 
band luminosity. Solid lines represent averages in five bins along 
the x-axis. The extent of the rms scatter is identified by dashed 
lines. The average error on individual points shown in the bottom 
right of each plot. 

light-ratio in old galaxies (solid lines) after correction for 
the mass-metallicity relation evident in Fig. [10] and Table 
[2] Similarly, the dashed lines in the right hand plots repre- 
sent the observed trends between luminosity and mass-to- 
light-ratios (solid lines) after correction for the luminosity- 
metallicity relation evident Table [2] The dashed lines can 
therefore be considered the relations for old galaxies of fixed 
metallicity. We note that the K band relations are almost 
unaffected by the corrections. This is due to the lack of sen- 
sitivity of K band luminosity to metallicity (see Table [TJ . 

The agreement between corrected log[M/L] values is ex- 
tremely good, with logarithmic slope against mass of ~0.15 
in all three bands. Against luminosity we find logarithmic 
slopes in the corrected values of -0.014, 0.010 and 0.031 in 
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Figure 12. Plot of mass-to-light-ratio in B, R and K bands 
against dynamical mass and K band luminosity. Thick solid lines 
reproduce the binned averages of the old galaxies from Fig. 1111 
Thick dashed lines represent these values corrected for metallic- 
ity. In B and R band plots, the K band profile is reproduced (after 
and offset ~0.9 dex) as a thin red line. The metallicity corrected 
values are remarkably consistent across the three wavebands. 



B, R and K bands respectively. With formal errors of ~0.01 
these are broadly consistent with a slope of zero. The re- 
markable consistency between the results for the three in- 
dividual wavebands suggest that the ages, metallicities and 
photometric predictions of the BC03 models are both reli- 
able and internally consistent, at least when they are used 
in the differential sense employed in this work. 

The results are also again in good agreement with Tru- 
jillo, Burkert & Bell (2004), who find slopes equivalent to 
0.20 and 0.14 in B and K bands respectively (compared to 
our 0.23 and 0.15). With the benefit of spectroscopy, we 
are able to clearly identify metallicity as the main contrib- 
utor to the stellar population effects. Our results are also in 
good general agreement with the previous studies of both 
Bernardi et al. (2003b) and Padmanabhan et al. (2004), 
who find similar trends in dynamical mass-to-light-ratios for 
large numbers of SDSS galaxies. 



2(1 - s 
a = — \ 

1 + 8 

and 



Where <I> is the average surface brightness within the ef- 
fective radius (R) and s is the observed slope in log[M/L] 
with log(Mass) from Table [2] The virial theorem for homol- 
ogous systems (i.e. in the absence of a tilt) would predict a 
= 2 and p = 1. 

The results of the comparisons between predicted a and 
and the (observational) literature values are shown in Ta- 
ble [3] The predictions using the observed slopes in mass- 
to-light-ratio with mass compare very well to the literature, 
values (consistent with the predictions to within ~2cr). We 
note that this is despite both our indirect approach as well 
as literature studies not having the luxury of isolating old 
galaxies. 

The comparisons indicate that, while most of the tilt 
of the FP is dynamical in origin, the mass-metallicity rela- 
tion in galaxies (Fig. I10p causes an increasing a and j3 with 
wavelengtlQ due to the decreasing sensitivity to metallicity 
of photometry at increasing wavelength. This trend culmi- 
nates in the observed FP tilt in the K band being almost 
identical to the underlying dynamical tilt in the FP (see pre- 
vious section). This is again due to the lack of sensitivity of 
K band luminosity to metallicity in old stellar populations, 
which makes the K band an excellent tool for probing mass- 
to-light-ratios in old (or uniform age) stellar populations. 

We have therefore shown our results to be consistent 
with a broad range of previous studies of galaxy mass- 
to-light-ratios and the Fundamental Plane. However, the 
conclusions we draw are not. To understand why this is 
the case, recall that the observed 'tilt' in the Fundamental 
Plane is often interpreted as a trend in mass-to-light-ratio 
with mass or luminosity. Indeed, a trend with mass of the 
form [M/L]txM s was the starting point of the previous 
consideration of the Fundamental Plane. However, careful 
consideration of the data show there is a flaw in this 
interpretation. The relation [M/L]cxM a can be simply 
rearranged into the form [M/LJrxlA This gives b=^^— . 
Therefore, for a=0.15, we would expect b=0.18, in stark 
contrast with the b=0 found. We must therefore conclude 
that the simple picture of a trend in dynamical mass-to- 
light-ratios with either mass or luminosity is excluded by 
our data. 



4.7 The fundamental plane 

While our data are not suitable for a direct analysis of 
the Fundamental Plane it is nonetheless interesting to com- 
pare our results with such studies in the literature. To this 
end we examined the implications of the varying slopes 
in mass-to-light-ratio with mass with waveband on fits 
to the Fundamental Plane (FP). The implications can be 
investigated using the measured slopes in log[M/L] with 
mass for old galaxies (Table [2] and Figs [11] and 11 2 p . Us- 
ing these values it is possible to predict the FP (given by 
log(R)=Qlog(cr)-|-/31og<I>-|-7) that one would expect to find 
using; 



In summary, we have shown that our observations are 
generally in good agreement with previous literature stud- 
ies. Furthermore, we have shown that the mass-metallicity 
relation with logarithmic slope 0.25 found in the data is 
sufficient to explain the differences in the observed trends 
in mass-to-light-ratios with mass and luminosity in B, R 
and K bands as well as the variations in the apparent 'tilt' 
found in previous studies of the Fundamental Plane. How- 
ever, once this stellar population effect has been removed, a 
significant (and consistent) trend in mass-to-light-ratio with 

2 Similar trends in a and (3 with wavelength were found by 
Bernardi et al. (2004) for the galaxy colours in the SDSS. 
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Figure 13. Mass-to-light-ratio in B, R and K bands against mass 
for both young and old galaxies. Average errors on data points are 
shown in each plot. Mass-binned averages of old and young galax- 
ies are shown as red and blue lines respectively. Blue dashed lines 
represent the mass-to-light-ratios expected in the young galaxies 
based on their central spectroscopic ages and mctallicities. While 
young galaxies deviate from the old population in a sense con- 
sistent with their younger ages, none have mass- to-light-ratios as 
low as their central ages would imply. 



mass remains in all three wavebands. Finally, we find that 
the data rules out simple models of variations in dynamical 
mass-to-light-ratio with mass or luminosity. Further analy- 
sis will require a detailed consideration of the distribution 
of galaxies in the full 3-dimensional space of the Fundamen- 
tal Plane. This is planned for a future paper. In this work 
we proceed by exploring the effects on photometry of young 
central spectroscopic ages. 



4.8 The effects of age on mass-to-light-ratio 

We now turn our attention to the effects of the young 
central populations detected in our 6dFGS spectra on the 
global properties of galaxies as indicated by their photome- 
try. There are 400 'young' (<1.5 Gyr) galaxies (see Section 
I4.5|l included in this analysis. 

Plots of mass-to-light-ratio with mass for both young 
and old galaxies are shown in Fig. 1131 Mass- binned averages 
of young and old data are shown as blue and red solid lines 
respectively. Also shown in these plots as dashed blue lines 
are model predictions of the effects of the differing ages and 
metallicities of the young and old populations in each mass 
bin. These are based on the comparison of luminosity values 
from the BC03 SSP models for the differences in age and 
metallicity of the young and old galaxies in each mass bin. 
A direct comparison of the data for young galaxies to these 
lines would therefore implicitly assume that the young cen- 
tral populations detected in the spectroscopy pervade the 



galaxies' populations as a whole - i.e represents the same 
fraction of the galaxy population at all radijfj. 

Young galaxies can be seen to be displaced with respect 
to the old galaxies in a sense consistent with their young 
ages. However, the displacements are small compared to the 
predicted values. They also become smaller towards higher 
masses. The data therefore show that, on average: (i) young 
populations do not pervade their host galaxies, but are in- 
stead centrally concentrated; ii) young populations must 
constitute relatively small fractions of total galaxy masses; 
and (iii) the mass fraction in the form of a young stellar 
population must decrease with increasing galaxy mass. 

To quantify these conclusions, for each mass bin, the 
fraction of total galaxy stellar mass and luminosity involved 
in the recent star-burst (fjvf and ft) were estimated. These 
estimates were made by a differential comparison of the 
BC03 predictions for mass-to-light-ratios to the observed 
values. In each waveband, the difference between the ob- 
served mass-to-light-ratios of the young galaxies from those 
of the old galaxies in the same mass bin (Alog[M/L] &s) 
was compared to the difference between model predictions 
for young and old populations of appropriate metallicity 
(Alog[M/L] mo£ j e ;). Denoting these mass-to-light-ratio differ- 
ences as s and r respectively, we derive: 



ft, 



10 -Alog[M/L] obs _ 1 1Q - 

or 

lQ-Alog[M/L] modcl _ I IQ- 



1 



(6) 



The derivation of Equation 6 is given in Appendix A. 
We estimate the luminosity fraction in the recent burst by: 



f L = 



rfiv 



(1 - f M ) + rh 



(7) 



Finally, the total mass (Mb urst , in solar masses) is esti- 
mated by: 



My 



Im.M„ 



(8) 



Where M 9a ; is the total mass of the galaxy. The deriva- 
tion of these expressions again assumes that the young pop- 
ulations at any given mass are seen against a background old 
population with a metallicity and mass-to-light-ratio appro- 
priate to that mass (as given by the mass-binned values for 
old galaxies in Table [2}. The young stellar population is 
nonetheless assumed to dominate the central regions (such 
that issues related to degeneracies with burst-strength do 
not arise). 

The estimates of mass- and luminosity-fractions and 
burst-masses are plotted against galaxy mass in Fig. 1141 
Despite the (expected) variation in luminosity fractions be- 
tween wavebands, agreement between the values for the 
mass fractions derived from the three wavebands is ex- 
tremely good. This suggests that both our assumption of the 
same underlying dynamical 'tilt' in young galaxies as exhib- 
ited by old galaxies is correct, and that our age and metallic- 
ity corrections are accurate. The trend shows a decreasing 
mass fraction with increasing mass falling from ~10% at 
10 10 ' 5 M s to 2% at 10 11 ' 5 M Q . This is in agreement with 
the study of Treu et al. (2005), who find the stellar fraction 



3 The analysis also implicitly assumes that the dynamical trends 
observed in old galaxies, whatever their cause, are also present in 



the young galaxies. 
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Predictions 


Literature 




Band 


a 




a 




References 


B 


1.25 


-0.81 


1.20 ± 0.06 


-0.83 ± 0.02 


j0rgensen et al. (1996) 








1.33 ± 0.05 


-0.83 ± 0.03 


Dressier et al. (1987) 


R 


1.29 


-0.82 


1.22 ± 0.09 


-0.84 ± 0.03 


Colless et al. (2001a) 








1.37 ± 0.04 


-0.825 ± 0.01 


Gibbons et al. (2001) 








1.38 ± 0.04 


-0.82 ± 0.03 


Hudson et al. (1997) 


K 


1.46 


-0.86 


1.53 ± 0.08 


-0.79 ± 0.03 


Pahrc et al. (1998) 



Table 3. Fundamental Plane coefficients. Values predicted from the slopes of mass-to-light with mass from Table [2] (see Section 14.611 
are shown in the left-hand column. These are to be compared with the values from the literature and values derived in this work. Good 
agreement is found. 
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Figure 14. The average luminosity (bottom) and mass (middle) 
fractions in recent bursts are plotted against average galaxy mass 
for the five mass bins of Figure [TT] Values derived from the B, R 
and K bands are shown in blue, red and black lines respectively. 
Also shown (top) are estimates of the actual masses (in Mq) 
involved in the recent bursts. A fixed 10% fraction is indicated by 
the dashed line. Agreement in mass estimates between wavebands 
is extremely good, suggesting an upper limit to the size of star 
bursts in these galaxies. 



and young (black points) central populations are shown. Av- 
erages (binned by luminosity) of the old and young pop- 
ulations are again shown as solid red and blue lines re- 
spectively. Slopes in the old galaxies are in good accord 
with literature values (thin black lines in each plot). As the 
sample includes late- type galaxies, particular amongst the 
young galaxies, these must lie close to the 'red sequence' in 
our colour-magnitude diagrams. This suggests they there- 
fore correspond to the 'dusty red-sequence' of Wolf, Gray & 
Meisenheimer (2005) . The lack of galaxies in the 'blue-cloud' 
is almost certainly the result of a number of selection effects. 
These include the effects of template mismatches during ve- 
locity dispersion measurements and the removal of galaxies 
with poor fits or large errors during age/metallicity mea- 
surement. 

Fig.[T5]shows that, for both colours, young galaxies tend 
to be bluer than the old galaxies. However, the offsets are 
extremely small. Our results from the previous section indi- 
cate that this is largely the result of the low mass fractions 
formed in the recent star formation events. Fig. [TS] there- 
fore shows that galaxies lieing on, or close to, the red se- 
quence may not be 'red and dead', but may still be forming 
modest numbers of new stars. Indeed, we note that most 
of the galaxies with young central populations also exhibit 
on-going star formation (Fig. [6}. 



formed in recent times varies from 20%-40% below 10 11 Mq, 
to below 1% above 10 115 Mq. 

Also shown in Fig. [14] are the estimates of average total 
mass involved in the recent burst of star formation. Agree- 
ment between wavebands is again good. However, the ab- 
solute values given here must be treated with caution, as 
they are calculated by the only non-differential application 
of BC03 models used in this work. Nevertheless, in a relative 
sense, the data exhibit a narrow range of burst-masses and 
an apparent upper limit on their size in high mass galaxies. 

4.9 Colour-magnitude relations 

As a final consideration we next investigate the implica- 
tions of our findings for colour-magnitude relations. Note 
that these relations are independent of dynamical effects. 

In Fig. [15] we plot the colour-magnitude relations B-K 
and R-K against M;(. Galaxies with both old (grey points) 



5 CONCLUSIONS 

We have measured ages and metallicities from the 6dFGS 
spectra of ~7000 galaxies of all morphological and emission 
types. We have demonstrated that the age-met allicity de- 
generacy has been broken in this study, particularly with 
regard to the very young and very old galaxies used in the 
subsequent analysis. These data for the central regions are 
compared with the global photometry from 2MASS and Su- 
perCOSMOS for the ~6000 galaxies for which such photom- 
etry is available. Combining spectroscopy and photometry 
we are able to investigate trends in global B, R and K band 
mass-to-light-ratios with the central stellar populations of 
the galaxies. 

An age-metallicity trend is identified in passive galax- 
ies that is consistent with the results of similar studies in 
the literature. The stellar populations of galaxies exhibit- 
ing emission lines (arising from either star-bursts or AGN) 
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Figure 15. The colour-magnitude relations B-K and R— K 
against M-k are shown for both young and old galaxies. 
Luminosity-binned average values are shown as blue and red lines 
respectively. Young galaxies are, as expected, bluer than old. 



follow a similar trend, but exhibit younger ages than their 
passive counterparts. 

To minimise the effects of complex mixtures of stellar 
populations, we confine our analysis to the data for very 
old or very young galaxies. We find a steep mass-metallicity 
relation in old galaxies (logarithmic slope of ~0.25), while 
in young galaxies the trend is consistent with zero slope. 

Using only the old galaxies, we find a dynamical trend 
in mass-to-light-ratios with mass of logarithmic slope ~0.15. 
This is in good agreement with values reported in the liter- 
ature. However, there is no agreement in the literature as to 
the cause of this tilt. Padmanabhan et al. (2004) concluded 
that most of the tilt is due to variations in dark matter 
content. From the modelling perspective, Robertson et al. 
(2006) perform simulations of a series of mergers and inves- 
tigate the scaling relations of their merger remnants. From 
mergers of gas-rich disk galaxies they derive a FP with a 
= 1.55 and j3 — 0.82, in good agreement with our K band 
values (see Table 3). In their models the tilt away from the 
virial scaling relation is caused by gas dissipative effects and 
not dark matter content. Trujillo et al. (2004) on the other 
hand attribute the tilt to 'non-homology', deriving almost 
identical contributions to the tilt from dynamical and stellar 
population effects as those we find in this work. Our find- 
ings, however, also show that whatever the cause of the tilt, 
it is not a simple in mass-to-light-ratio (or its scatter) with 
mass or luminosity. Instead we find a that a more complex 
description will be required; this will be a subject of a future 
paper. 

Considering extremely young galaxies, our results show 
that the properties of the central stellar populations iden- 
tified by the spectroscopy do not generally correlate with 
the global properties of their host galaxies. Indeed, the data 
indicate that recent bursts of star formation are limited to a 
mass of <10 10 Mq and are strongly centrally concentrated. 



Consequently, in high mass galaxies exhibiting young cen- 
tral populations, recent star forming events have negligible 
effect on their global photometry. As a consequence we have 
shown that galaxies sitting on, or near, the red sequence 
may not be 'red and dead', but may still be forming modest 
numbers of new stars. 

We therefore conclude that the young ages found in 
studies of galaxy centres in this and other spectroscopic 
studies should not be considered representative of the age of 
their host galaxies as a whole. They may, however, be reason- 
ably interpreted as indicators of recent assembly events. As 
a result of these considerations, the small scatter in photo- 
metric scaling relations, which initially motivated the mono- 
lithic collapse models, ceases to be inconsistent with either 
the hierarchical merging model or the results of the Lick 
index studies indicating young central populations. 

In future work we will include the results from the 
completed 6dF Galaxy Survey and report on the full 
chemical properties (including a-abundance ratios) of the 
central regions of galaxies. The effective radii, peculiar 
velocities and environmental parameters of the 6dFGS 
sample, currently being analysed by the 6dFGS team, 
will be included. This will allow not only more accurate 
determination of photometric properties and a more de- 
tailed analysis of the Fundamental Plane, but also the 
investigation of the strong quantitative predictions of trends 
in galaxy formation/assembly epochs and metals-content 
with galaxy mass and environment. 
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1 Appendix 1. Mass fraction derivation 

The logarithm of the mass-to-light ratio in a galaxy compris- 
ing two populations; one young and one old, can be written: 

log[M/L] obs = log[- ]. (1) 

-L'old ~r J-Jyoung 
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This parameterisation assumes that the mass of the 
galaxy is unaffected by the starburst that generates the 
young stellar population (the galaxy converts some of its 
own mass into stars). 

Now, expressing the fraction of the observed stellar 
mass in young stars as fjf and the ratio of mass-to-light- 
ratios of young and old populations as r, we can write: 

log[M/L] obs = log[ M ]. (2) 

(1 — IM JJ-Jold + riM-Lold 

Rearranging: 

log[M/L] obs = log[M/L] old - log((l - f M ) + rf M ). (3) 
Denning: 

A log[M/L] obs = log[M/L] obs - log[M/L]oid. (4) 
We get; 

A log[M/L] obs = -log((r - l)f M + 1). (5) 
Rearranging: 

1Q -A log[M/L] ob3 = (r _ 1)fM + ^ (6) 

and; 

1Q -A log[M/L] ob3 _ ^ 



f M 



(7) 



The ratio of mass-to-light-ratios (r), as denned above can 
be estimated using the BC03 models of appropriate age and 
metallicity: 

r _ l^oglM/LlBCOSold-lloglM/LlBCOSyoung 



I.e: 

r = 10 



-A log[M/L], 



(8) 



(9) 



We thus derive the fraction of the stellar mass in the 1 Gyr 
old population as: 

10 -Alog[M/L] ob 



fM = 



1 



10-Alog[M/L] modcl _ I 



(10) 



We note that, by the differential use of both the 
observations and models, mass fractions are derived without 
knowledge of the true stellar to dynamical mass ratios. 
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